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Abstract 

The potential energy surface (PES) describing the interactions between Li2( 1 S+) and 4 He and an exten- 
sive study of the energies and structures of a set of small clusters, Li2(He) n , have been presented by us in a 
previous series of publications fl, 2, 3]. In the present work we want to extend the same analysis to the case 
of the excited Li2(a 3 £+) and of the ionized Li^ (X 2 £+) moiety. We thus show here calculated interaction 
potentials for the two title systems and the corresponding fitting of the computed points. For both surfaces 
the MP4 method with cc-pV5Z basis sets has been used to generate an extensive range of radial/angular 
coordinates of the two dimensional PES's which describe rigid rotor molecular dopants interacting with one 
He partner. 
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I. INTRODUCTION 



The accurate evaluation of the interaction potential between a dopant molecule or ion and a 
single He atom is an important step in the study of small-size doped helium clusters because it 
can be used to set up the total potential acting within them, an essential first step in determining 
their structure and dynamics. Nano-sized helium droplets in which a molecular or an atomic 
impurity has been added by means of picking up techniques provide, indeed, a unique environment 
for high precision spectroscopic studies of the molecular solute [4] especially for weakly bound 



molecular species O, la LZD that are otherwise difficult to analyze. Among the many dopants that 
have been experimentally and theoretically analyzed, the alkali metal atoms and dimers have many 




interesting properties: from experimental O, LZl IBO and theoretical HI 121 12D evidence it turns out 
that they normally reside on the surface of an 4 He droplet, usually forming a slight "dimple" on 
that surface. When two or more alkali atoms are attached to the droplet they eventually meet on its 
surface forming a singlet or a triplet molecule. Spectroscopic identification of the molecules shows 
lO] that dimers in the triplet state outnumber the singlet molecules by a factor up to 10 5 . This is 
probably due to the differences in stability between the triplet dimers and the singlet molecules 
Bo] whereby the energy release due to the formation of a singlet molecule is more likely to lead to 
the detachment of the molecule and to the almost complete disappearance of this kind of dopant. 
A similar behaviour, i.e. one in which high spin states are preferred, has been noticed also in 
the formation of small clusters of alkali atoms grown on helium droplet surfaces ifiol : the very 
small binding energies to the surface allows for a spontaneous "selection" of the cluster in higher 
spin states (it is worth pointing out that, in these experiments, lithium atoms were shown to form 
only dimers while heavier alkalis produced also larger structures). Very recently experiments have 
been also performed on alkali atoms in 3 He droplets Jill where both theoretical predictions and 
experiments [ 12] have shown a very similar behavior to that of the bosonic moieties. The possible 
formation of cold heteronuclear molecules has attracted the attention of researchers in the last few 
years (see for example the special volume edition in which Ref. Jl^l can be found) because of 
the possibility of studying the effects of polar interactions in cold or even quantum degenerate 
molecular samples: helium droplets may therefore allow the spectroscopic study of weakly bound 
(triplet) heteronuclear dimers. 

In all the experiments mentioned above, the helium droplet is used as a convenient, but inert 
matrix in which weakly bound molecular species are accumulated and studied. Helium droplets, 
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however, can also be used as microscopic cryo-reactors which actively initiate series of chemical 
reactions. Very recently, in fact, it has been shown possible to follow chemical reactions inside 
a helium droplet using, among others, simple diatomic dopants. The chain of reactions is 
usually initiated by the ionization of one of the helium atoms [14] with an ensuing charge transfer 
I15I flfii [nl process that moves the charge onto the impurity because its ionization potential is 
usually lower than He. The impurity gives rise, in turn, to a variety of further, secondary chemi- 
cal reactions: the excess energy due to any exothermic reaction is rapidly thermalized within the 
droplet itself by helium atom evaporation. The rapid quenching of the internal degrees of freedom 
of the partner species afforded by the superfluid environment of the droplets has a profound effect 
on the various branching reactions and, in some cases, can ultimately stabilize different interme- 
diate complexes which would instead only remain as transient species in the gas phase. The latter 
condition therefore becomes an interesting feature of such nanoscopic cryostats in the sense that 
any further step initiated by ionization of the cluster can follow specific pathways which are likely 
to be different from those under the corresponding gas phase conditions. 

Concerning the experiments involving neutral dopants, we believe that the microscopic struc- 
ture of the helium surface-dimer complex is a key piece of information for the full understanding 
of the experimental data. In this work we continue the analysis of the lithium-He system by pre- 
senting an accurate Potential Energy Surface (PES) for the dimer in the triplet state Li 2 ( 3 E+) and 
one Helium atom. 

When moving to the problem of ionic dopants, and therefore to the fascinating possibility of 
conducting ionic reactions at low temperatures in the droplet, the study of the microscopic struc- 
ture of the solvation site around the ionic molecule may also help the understanding of the possible 
reactive pathways which are eventually followed. This is the reason why, together with the pre- 
viously mentioned PES, we also want to present here similar calculations on the Li^ ( 2 £g ) — He 
system. These two PES's may be considered as the building blocks of the complex interactions at 
play in the He droplets or in smaller helium clusters: they describe two rather extreme situations of 
either a weak, dispersion-type interaction or of a much stronger, markedly orientational ionic PES. 
Although the forces at play in the large systems may be extremely difficult to compute, a simpli- 
fied approach where the total potential is obtained as a sum of two-body potentials represents a 
realistic route because the three-body forces and the higher order terms, especially for the neutral 
system, provide fairly negligible contributions to the total energies. For their possible relevance in 

no 

ionic systems see e.g. Refs. [18] and [19] for further details. 
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II. THE AB-INITIO CALCULATIONS 



Both surfaces have been calculated using the MP4(SDTQ) method (without freezing the lithium 
core) with the cc-pV5Z basis set employed within Gaussian03 All energies have been cor- 
rected for BSSE using the Counterpoise procedure I21L Since, in either case, the distortion due to 
the helium atom over the electronic structure of the molecule is rather weak, we have decided to 
keep the molecule at its equilibrium geometry. The optimal geometry of the molecular species was 
determined by a separate MP4 optimization on both the triplet Li 2 and ionic Li^ . The internuclear 
distances that we have used are r = 4.175 A and r = 3.11 A for Li 2 and Li^ respectively. For 
both surfaces we have used a grid of Jacobi coordinates with different radial geometries and 19 
angles from to 90. For the Li^ — He case we have calculated for each angle a number of points 
between 50 and 64 (2.5 < R < 12.0 A ) for a total number of 1,135. For the Li 2 — He system the 
number of radial points for each angle varied from 30 to 80 (2.0 < R < 14.0 A ) for a total of 847 
points. 



III. THE FITTING PROCEDURE 

In order to make it easier to employ the two potential energy surfaces in further studies,we 
have decided to obtain a suitable analytical fitting of the raw points using a non-linear fitting 
procedure based on the minimization of the square deviation, as obtained by means of the efficient 
Levenberg-Marquadt method J3]. 

The full interaction can be written as 

v tot = v(R a , e a ) + v(R b , e b ) + v lr (r, e) (i) 

where the coordinates are those reported for clarity in Figure 1. The first two contributions repre- 
sent the anisotropic interactions at short range and are written in terms of Legendre polynomials 

nmax Imax 

V(R a , 0a)=J2Yl R - e M-PRa) ■ Pl{cOsd a ) ■ C a nl (2) 
n=0 1=0 

where, given the symmetry of the homonuclear molecule = C h nl . The long-range contribution 
is instead expressed in Jacobi coordinates and given by 

V LR (R, 6) = £ £ ■ P ^ cos ^ ■ « <3) 



JV L=0 




FIG. 1 : Coordinates used in the fitting formula 
where the fjy damping functions are those defined within the well-known Tang-Toennies empirical 



potential modeling [23]. 

For the neutral system (Li 2 -He) we have employed nmax = 4 and Imax = 6 and we have 
limited the long range expression to the value N = 6, i.e. the latter was simply given by two terms 

V LR = f 6 ((3R) ■ R- 6 ■ [Po(cos8)C™ + P 2 (cos6)C^ R }. 

This very simple analytical representation of the entire set of unweighted 847 ab-initio points 
yielded a standard deviation of 0.0046 cm -1 . For the ionic system (Li^-He) we have employed a 
more flexible representation given by nmax = 8 and Imax = 4 and two long range anisotropic 
terms: C4 and Cq. The latter terms were employed in order to represent correctly the long range 
multipolar expansion of a neutral, polarizable atom interacting with a point-like charge ion. The 
long range expansion therefore becomes: 

V LR = U(J3R) ■ R~ A ■ [Pv(cos6)C^ R ] + f 6 (PR) ■ R~ 6 ■ [P Q {cos6)C^ R + P 2 {cos6)C^}. 

The total number of 1,135 points was reduced to 1,065 by excluding the repulsive energies with 
values above 5,000 cm -1 . A simple weight function W(E) was used in order to obtain a better 
representation of the interaction region and to reduce the importance of the repulsive geometries. 
Hence we drop the weight of the points in the less important highly repulsive regions of the inter- 
action 

[l for£< 1,000 

W(E) 

IE]- 1 for 1,000 < E < 5,000 
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FIG. 2: Potential energy curves for Li2( 3 £)-He at 9 = (solid line) and = 90(dashed line) as obtained 
from our fitting. Filled-in circles are ab-initio energies. 

The final standard deviation was 0.073 cm -1 . As an example of the accuracy of the two fittings 
we report in Figure 2 (neutral system) and Figure 3 (ionic system) two selected cuts through the 
PES's at 9=0 and 9=90 and compare them with the ab-initio points. As is evident from the figures, 
we are able to describe with good accuracy both the attractive region and the lowest section of 
the repulsive wall. In either case, however, we were not able to provide reliable fitting values for 
R < 1.0 A . All the coefficients are reported in Table I and Table II. FortranW subroutines are 
available on request from the authors. 




our fitting. Filled-in circles are ab-initio energies. 
IV. RESULTS AND DISCUSSION 
A. The Potential energy surfaces 

A 3D representation of the two PES's is presented in Figure 4 as isoenergetic contour plots. As 
can be seen from the figure, the two PES's are completely different in shape and in strength, as 
one could easily have expected. The triplet interaction is very weak (very similar to the one we 
have already calculated for the singlet molecule in Ref. tl|,|3j]) and has its absolute minimum at 
9=90 with an interaction energy of ~ —2.4 cm -1 . This is a very different anisotropy with respect 
to that exhibited by the singlet dimer JilQ], where the collinear configuration provided the global 
PES minimum energy and it also exhibits deeper well values than those of the former case. The 
extreme weakness of the interaction of Li 2 with He is indirectly supported by the experimental 



7 




FIG. 4: Contour plots of the PES's for Li 2 ( 3 S+)-He (upper panel) and Li^( 2 S+) - He (lower panel). 
Dashed contours are positive energy isolines. 
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TABLE I: Coefficients for the Fitting of Li 2 ( 3 S)-He 



n i 




n 1 


/- 1 


n 1 




n 1 




00 


22763.19 


1 3 


-122598.06 


26 


-825.49 


42 


1291.90 


1 


167466.30 


1 4 


57319.62 


30 


-2718.38 


43 


-513.93 


02 


-171109.44 


1 5 


-11157.70 


3 1 


1618.93 


44 


69.19 


03 


116907.07 


1 6 


567.34 


32 


672.89 


45 


16.31 


04 


-44321.81 


20 


-5532.20 


3 3 


-2363.22 


46 


-7.15 


05 


-675.99 


2 1 


31713.40 


34 


2099.82 


NL 


NL 


06 


3002.69 


2 2 


-43241.69 


3 5 


-816.34 


60 


-560011.07 


1 


5967.45 


2 3 


37506.82 


3 6 


157.05 


6 2 


105854.86 


1 1 


-143899.90 


24 


-20326.82 


40 


648.35 







1 2 


162737.35 


25 


5713.51 


4 1 


-1524.99 




1.09998326 



findings [7] which indicate that triplet molecules also reside on the surface of the droplet and are 
not efficiently solvated by the helium atoms. 

Given the weakness of the calculated interaction it would be also be interesting to provide a 
modification of it that makes use of semi-empirical guesses for the two-body contribution due to 
the Li-He interactions. We can, in fact, modify the calculated interaction energies by employing 
2-body semiempirical potentials along the lines which we have already followed for the singlet 
molecule in Ref. |3]. It is, in fact, possible to produce rather simply a semi-empirical potential 
V'(R, 6) by using the formula 



V(R, 



V(R, 6) - V ab (R a ) - V ab (R b ) + V sm (R a ) + V sm (R b 



where V sm (R a ^) is the semi-empirical potential proposed by Toennies and coworkers 12311 for the 
dimer Li-He, V ab (R a> b) is a simple fitting of the ab-initio points calculated by us at the same level 
of accuracy (MP4/cc-pV5Z) for the Li-He diatomic curve and R a , R b are the two Li-He distances 
within the trimer. Our ab-initio data, the fitting curve and the model potential of Ref. B23I1 for the 
Li-He pair are reported in Figure 5. The resulting full potential energy surface V'(R, 6) calculated 
with the model potential is stronger than the one from purely ab-initio data, although the increase 
in well depth is still not sufficient to allow for possible solvation of the molecular impurity inside 



TABLE II: Coefficients for the Fitting of Li 2 ( 3 S)-He 



n 1 




n 1 




n 1 




n 1 




00 


1794639.56 


23 


-26905991.00 


5 1 


7206051.75 


74 


-4773.84 


1 


-4891313.47 


24 


-5708309.41 


5 2 


-7598549.02 


80 


4440.63 


02 


9692348.34 


30 


-7120382.44 


5 3 


4029614.92 


8 1 


-12773.88 


03 


-6003108.42 


3 1 


33556992.50 


54 


440943.65 


82 


12647.36 


04 


-991282.46 


3 2 


-41302351.90 


60 


529353.73 


8 3 


-7220.96 


1 


-3645612.91 


3 3 


21971925.50 


6 1 


-1650913.80 


84 


760.89 


1 1 


20414626.80 


3 4 


4504104.73 


6 2 


1675760.54 


NL 


^NL 


1 2 


-33068152.00 


40 


5071552.76 


6 3 


-909621.15 


40 


-10627.48 


1 3 


19171419.60 


4 1 


-19611375.40 


64 


-32113.64 


60 


+43087.50 


1 4 


3756985.86 


42 


22000474.60 


7 


-73683.34 


62 


-212414.50 


20 


6019745.96 


4 3 


-11557831.20 


7 1 


217607.26 




P 


2 1 


-35226804.20 


44 


-1959452.27 


7 2 


-216505.97 




3.01084652 


22 


48899465.50 


50 


-2124762.77 


73 


120820.53 







n 

the helium clusters (see also our similar conclusions for the singlet Li 2 in Ref. OD). 

The ionic interaction reported in Figures 3 and 4, on the other hand, is much stronger and 
presents its minimum at 0=0. Here the situation is completely different and, although no calcu- 
lations have been carried out as yet for the cluster structures, it is not difficult to see that the Li^" 
impurity would be strongly solvated in liquid helium droplets and therefore it is likely to get local- 
ized at the droplet center with a solvation shell strongly bound to it (for analogous situations one 
should look at the existing results on the ionic Li or Na impurities described in Refs. JlslQ])- 

B. Estimating the 3-body effects in the ionic system 

As we have pointed out in the introduction, in order to treat small and medium sized clusters 
with 2-50 helium atoms, one should be able to rapidly evaluate the total potential energy acting 
within them. One of the most used and successful approaches consists in approximating the total 
interaction as a sum of 2-body terms, initially neglecting any non-separable 3-body contribution. 
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FIG. 5: Potential energy curves for Li-He pair. Filled in circles: present ab initio data; solid line: model 



potential of Ref. 12311 : dotted line: fitting of calculated ab-initio data using the following formula:y(i?) = 
En=o a n R n e' bR - h(bR)C 6 R~ 6 with these parameters: a =271476, ai=-33951.3, a 2 =-1.11528-10 6 , 
6=1.76229 and C 6 =77803.7 

While this has been proved to be a very accurate procedure for doped helium clusters with a neutral 
impurity, it may represent a source of error in ionic clusters where non-separable interactions 
among the induced multipoles in the first solvation shell may be important. However, it has been 
recently shown in various works on ionic dopants in rare-gas clusters including helium that these 
effects are small and should not, when included, alter substantially the geometries or the energies 
of the clusters (for a general discussion with anionic dopants see refs. 11911 and [25, 26], for 



positive ions see Refs. 118 



27]). 



In order to verify the applicability of the sum-of-potentials approximation we have carried out 
at a consistent level (cc-pV5Z/MP4) a series of ab-initio calculation on the Li^He 2 system for 
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which such 3-body effects may arise. The analysis carried out here is analogous to what we have 



already reported in Ref. [27] for the LiH + dopant. The 3-body forces should mainly originate from 
the induced multipoles (attractive and repulsive contributions) and by the weaker Axilrod-Teller 
effects tfioll . For our preliminary study, we have therefore chosen four different geometries of the 
complex where the two helium atoms are close enough to the ion molecule in order to contribute 
significantly to 3-body forces. The first geometry (geometry A) is weakly repulsive (~ 70 cm -1 ) 
with the two helium atoms at 3.11 A and is the upper one reported in the inset of Table 3. The 
second (B) and the third (C) geometries are similar in shape to the second one sketched in the same 
table, but they differ for the values of the distances between the atoms. The last one (geometry opt) 
is also similar in shape to the latter, but it comes instead from a full minimization at the MP2/cc- 
pv5Z level of the entire complex: all the relevant distances are reported in Table 3. The analysis of 
3-body forces is done in the following way: for each geometry we calculate the interaction energy 
of the entire Li^He 2 complex, of one of its 2-body fragments (there are two identical ones in 
symmetrical geometries) Li^He and of the remaining He 2 dimer by using the following formulae 
where the geometry is fixed, it is the same in each fragment and the counterpoise correction is 
used for each interaction energy: 

y[Li+He 2 ] = £[Li+He 2 ] - £[Li+] - 2E[He] 
F[Li+He] = £[Li+He] - E[Li+) - E[Re) 

V[He 2 ] = £[He 2 ] - 2E[Re] (4) 

the residual 3-body interaction is simply calculated by using the expression 

y[Li+He 2 ] - 2y[Li+He] - V[Ee 2 ]. (5) 

As can be seen from Table 3, the 3-body interaction is always a small percentage of the total 
interaction in the various geometries except for geometry A where it represents more than 10% of 
the total interaction. It is however important to note that the geometry A corresponds to a repulsive 
geometry for the Li^-He fragment and therefore it is not relevant for optimization purposes. 



V. CONCLUSIONS 



We have computed two accurate Potential Energy Surfaces for two different systems that are 
of interest for the experimental and theoretical study of helium droplets doped with alkali metal 
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TABLE III: Distances (A) and energies (cm 1 ) of the test geometries used to estimate the importance of 
3-body effects in the ionic system. Energies have been corrected for BSSE. 

^ Rillo He) ^ 
R(Li He) 




Quantity 


Geom A 


Geom B 


Geom C 


Geom opt 


R(Li-Li) 


3.11 


3.11 


3.11 


3.04382 


R(Li-He) 


3.00 


3.00 


3.60 


1.91530 


R(He-He) 


3.11 


3.00 


1.25 


2.52866 


F[Li+He 2 ] 


66.113 


-223.547 


11376.770 


-811.745 


^[Li+He] 


38.745 


-110.088 


-52.863 


-404.729 


^[He 2 ] 


-3.436 


-2.9471 


11367.744 


22.8860 


V(3B) 


-7.941 (12%) 


-0.424 (0.2%) 


114.752 (1%) 


-25.173 (3.1%) 



molecules. The two molecules considered here are the triplet state of Li 2 and the ground state 
ionic Li^ . As should be expected, we found the two PES's to be markedly different: 

• the interaction of the neutral moiety is similar to the one we have already studied for ground 

nn 

state singlet Li 2 111121]: a very weak interaction that confirms once more the tendency of high 
spin compounds of alkali metals to reside on the surface of helium droplets. A modified ver- 
sion of the same interaction that gives rise to slightly deeper potentials can be obtained fol- 
lowing our earlier proposal in Ref. [3] and also found to yield weaker interaction potentials 
than those between He partners. 

• the triplet dimer interaction with helium is showing here a markedly different anisotropy 
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from that found earlier on for the singlet state of Li 2 : the minimum energy configuration 
is, in fact, given by T-shaped structures as opposed to the linear structures obtained for the 
singlet interaction. 

• the ionic interaction, instead, is much stronger and more orientation-dependent: it should 
therefore lead to full solvation of the molecular moiety inside the droplets. 

• a preliminary analysis of the three-body effects on the interactions with more He atoms 
indicates that such effects are relatively small and should therefore allow the use of an 
approximate description of the full potential energy landscapes in Li^~He n clusters in terms 
of two-body potentials. 

Both surfaces have been fitted using rather simple analytical expression reported in the present 
paper and which therefore provide working quality interactions potentials for the title systems. We 
believe that such potentials are an important step in the modelling of the larger clusters behavior 
because they may be used to set up the total interaction for much larger systems whenever using the 
sum-of-potential approach is found to be a realistic alternative: our present study and preliminary 
analysis do seem to suggest that this may be the case for both the present system. We shall verify 
such a possibility via our ongoing calculations for the doped 4 He clusters. 
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